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Abstract 

By the methods of the time series analysis we investigate human vestibular 
and optovestibular nystagmus time series. They can be stationary or nonsta- 
tionary and are between the classes of purely periodic or purely chaotic time 
series. The amplitudes of the autocorrelation functions of the time series have 
an unexpected peak between 7 and 9 seconds on the time axis. The singular 
spectrum analysis shows that the important information about the time series 
is concentrated in the first four to six principal components. The influence 
of the sensory input modality, which can be changed by opening or closing 
eyes of the investigated persons, leads to changes in the histogram, but leaves 
almost unchanged other important characteristics of the time series such as 
the statistical dimension or the positions of the dominant frequencies in the 
power spectrum. Thus the behavior of the studied time series shows a pres- 
ence of nonlinearities which could play an important role for realization of the 
nystagmus phenomenon. 
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Nystagmus [^] is a rhythmic involuntary eyes movement that stabilizes gaze in the space 
during head movements or rotations of the environment. By means of this physiological 
mechanism the images on the retina do not move with respect to the main optic coordinates 
* corresponding author 



1 



and a clear vision is maintained. The nystagmus reactions are part of the general function 
with critical biological importance - orientation in the space. These reactions are very 
important for health organisms under the action of extreme influences such as the sport, 
astronautics, sailing, construction of high buildings, high-speed transport, etc. In all of 
these cases the professional selection based on medical examination include as a main 
factor the nystagmus reactivity. 

Many structures are involved in the realization of the nystagmus reaction: from 
vestibular and optovestibular receptor structures to nervous-oculomotor apparatus with 
involving of large brain zones (from brain-stem to cortex) in the sensory-motor integration. 
As a final product the nystagmus is a carrier of important information for the functional 
status of the forming structures. This makes the nystagmus an obligatory criterion in 
the clinical diagnostic of neurological, vestibular and eyes disfunctions. The nystagmus 
eyes movements are to the same amount but in the opposite direction with respect to the 
direction of the movement of the head. The nystagmus has two alternative phases of slow 
and fast movements in opposite directions. The slow movement phase is connected to the 
graceful following eye movements and contains the biological part of the nystagmus reflex. 
The fast movement phase is connected to the jump-like saccades and reflects only the me- 
chanical reposition of the eyes. The vestibular and the visual system ensure the afference 
of the nystagmus called vestibular or optovestibular in dependence on the sensory input 
modality. With respect to its external manifestation the two kinds of the nystagmus do 
not differ each from another. 

The goal of this paper is to investigate characteristic peraccelatory nystagmus time 
series by the methods of the nonlinear time series analysis [^,^,'^]. Despite the progress in 
the area of the registration and evaluation techniques all characteristic nystagmus quan- 
tities have large deviations from its assumed physiological normal vahics. Thus it is very 
difficult to distinguish between the norm and the pathology. This makes the problem for 
the investigating of the nystagmus highly actual today. The experimental data are chosen 
among the time series, obtained under physiological conditions from a healthy volunteer, 
a man aged 32 years, without history of past or present vestibular, otologic, opthalmo- 
logic or neurologic diseases, that potentially could affect the nystagmus results. For the 
experiment we have used an electronically controlled rotating chair (Toennies, Freiburg). 
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Stimuli consisted of acceleration steps (ib8°/s^, 190°/s) to the right and to the left. Ro- 
tation test have been done under two different visual conditions: in complete darkness 
with eyes closed (vestibular nystagmus), and with eyes open in a lighted optokinetic drum 
(optovestibular nystagmus). During passive body rotation in the light concomitant visual 
and vestibular information interact in order to provide better ocular stability [^]. Hori- 
zontal eye movements have been recorded bipolarly from surface electrodes at the outer 
canti by means of conventional electronystagmography. It is based on the registration of 
the biopotential difference which arises between the retina and the cornea in the cases 
of different movement of the eyes. The eye-ball is a dipole: the one pole, the cornea is 
positive charged, and the another pole, the retina is negative charged. When the human 
sees in the straight direction the potential between the cornea and the retina is small. 
When the sight changes its direction this potential can increase ten times. Prom the reg- 
istered nystagmus reaction we can estimate the basic nystagmus parameters as length in 
the time, frequency, amplitude and velocity for the different components of the nystagmus 
cycle as well as for the entire nystagmus reaction. Normally one calculates the velocity of 
the slow phase of the nystagmus. It is a functional analogy of the compensatory eye-ball 
movements. 

The nystagmus time series can be stationary or nonstationary ones. This is shown in 
panels a and b of Fig. 1. The time series have complicated structure because of the influ- 
ence of afference sense structures, effectors occumotor apparatus and connecting segments 
of the central nervous system. The stationary time series for vestibular nystagmus are 
shown in panel a. The time series connected to the optovestibular nystagmus shown in 
panel b are clearly nonstationary with a trend given by the relationship V = 8.3 — 0.045i, 
where V is the angular deviation of the eye-balls from the strait sight direction in the 
horizontal plane and t is the time. After a detrending these time series become stationary 
one and we shall use this stationary form in the analysis below. We perform time delay 
embedding of the time series, obtaining the value of the delay (r = 4) and the minimum 
embedding dimension (m = 4) by means of the first minimum of the mutual information 
P] and the method of the false nearest neighbors [^] . The panels c and d of Fig. 1 show 
the results. The obtained manifolds contain points scattered almost on one plane and we 
can expect that the correlation dimensions D2 ^] of the investigated time series are 
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between 1 and 2. Indeed D2 = 1.3 it 0.14 for the vestibular nystagmus time series and 
D2 = 1.4 lb 0.12 for the optovcstibular nystagmus time series. Thus the change of the 
influence of the environment do not change quite significantly the value of the correlation 
dimension of the time series and they remain in the middle class between purely periodic 
and purely chaotic motions in the phase space. 

Fig. 2 shows the autocorrelations, histograms and the lower parts of the power spectra 
for the investigated time series. The autocorrelation function exhibits a periodic behavior 
with large region of almost exponential decaying of the extrema of A{t). After this the 
correlation of the data increases, and the beginning of the regions of increasing are almost 
on the same place in the two time series. The autocorrelations of periodic processes 
are periodic and the autocorrelations for chaotic processes decay exponentially . The 
autocorrelations of the investigated time series show exponential decaying of the amplitude 
of the maxima for small r and thus are close to the chaotic systems. But between 350 
and 450 sampling units we observe increasing of amplitude of the autocorrelations for 
both time series. This peak between 7 and 9 seconds of time hints that the nystagmus 
dynamics is nonlinear one. For stationary time series the power spectrum does not change 
in course of evolution of the time series. The power spectrum does not show the famous 
1 //- behavior. The histograms shown in panels c andrf of Fig. 2 show that the distribution 
of the values of the time scries is far from the Gaussian one. The vestibular nystagmus 
data exhibit fat tails and the optovcstibular data histogram has two peaks. The power 
spectral density is concentrated in the same two frequencies for the vestibular and for the 
optovcstibular nystagmus. Thus the changing environmental conditions do not change the 
dominant frequencies of the nystagmus dynamics. 

The singular spectrum and the principal components of the time series are convenient 
tool for analysis of stationary and nonstationary time scries . ^^]. The largest singular 
values of the investigated time series are connected to their periodic components. The 
contribution of each principal component to the total variance of the time series is a 
measure of the important information about the series which is stored in the correspondent 
principal component. Thus by means of the variance we can determine the number of the 
important principal components and we are able to construct on their basis noise-free 
time series, which have all important features of the original series. The number of the 
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important singular values gives the statistical dimension of the time series which is an 
upper bound on the correlation dimension. The results of the singular spectrum analysis 
performed in twenty-dimensional phase space of vectors, obtained from the time series by 
means of unit delay, are presented in Fig. 3. Panels a and b show the singular values. 
Usually the nonlinear oscillations are presented by couple of nearly equal singular values. 
We can recognize three such oscillations in panels a and b of Fig. 3. Another observation 
is that the first four singular values are much larger than the other ones. Thus we can 
conclude that the statistical dimension of the time series is Sy = Sg = 4. The panels 
c and d present the percents of the total variance of the time series contributed by the 
corresponding principal component. As we can see almost all important information in 
the investigated time series for the vestibular and for the optovestibular nystagmus is 
concentrated in the first four principal components, i.e. the change in the sensory input 
do not change the number of important principal components. We can project the time 
series into the subspace of the eigenvectors corresponding to the largest four to six singular 
values. As a result we shall have time series which contain all important information of 
the original time series and are free of the influence of low-amplitude and high-frequency 
components such as noise. The remaining panels of Fig. 3 present the first five principal 
components for the time series for the vestibular nystagmus. 

As concluding remark we note that the obtained in this article vestibular and op- 
tovestibular nystagmus time series belong to the class of time series which are between 
the classes of purely periodic and purely chaotic motions. The interests in investigation of 
such time series shows a steady increasing in the last years. The autocorrelations of the 
time series show sign of both periodic and chaotic behavior and in addition they present 
evidence for nonlinearity of the time series by the existence of large local peaks of the 
amplitude after region of exponential decreasing of the peak amplitudes. The form of 
the histograms of the time series differs from the Gaussian one and the most important 
frequencies in the power spectrum are the same for both kinds of investigated nystagmus 
time series. The singular spectrum analysis shows that the effects of low-amplitude and 
high-frequency components of the time series can be easily eliminated by projecting the 
series onto the subspace of four to six appropriate eigenvectors. Finally, if stationary time 
series are result of purely linear process, they can have periodic, exponentially increasing 
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or exponentially decreasing behavior in the case of lack of noise. Other kind of behavior 
are evidence for presence of noise or nonlinearity. Noise do not affect significantly the 
investigated time series - several schemes for noise reduction did not changed significantly 
the properties of the time series (as it is known the nonlinear characteristics are very 
sensitive to presence of noise). The closeness of the values of the various characteristic 
quantities when the eyes are open or closed shows that the nystagmus time series are 
robust against changes in the sensory input modality. 
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FIGURE CAPTIONS 



1. Fig.l 

Time series and their images in the reconstructed phase space. Panel a: vestibular 
nystagmus time series. Panel h: optovestibular nystagmus time series. One unit for t 
is 0.02 s. One unit for V is 0.6 degrees deviation of the eyes in the horizontal plane for 

panel a and 0.9 degrees deviation for panel h. Panel c: image of the time series in the 
subspacc of the first three coordinates X, Y, Z of the four-dimensional reconstructed 
space of the delay vectors for the case of the vestibular nystagmus time scries. Panel 
d: image of the time scries in the subspace of the first three coordinates X.Y, Z of 
the four-dimensional reconstructed space for the optovestibular time series. 

2. Fig.2 

Autocorrelations A, histograms H and power spectra P for the nystagmus time 
series. Panels a,c,e: quantities for the vestibular nystagmus time series. Panels 
b, d, e: quantities for the optovestibular nystagmus time series. One unit for A is 25 
Hz. 

3. Fig.3 

Results from the singular value decomposition and principal component analysis of 
the time series. Panel a: singular values for the vestibular nystagmus time series. 
Panel h: singular values for the optovestibular nystagmus time series. Panel c: Por- 
tion of the total variance, due to the corresponding principal components for the 
vestibular nystagmus. Panel d: Portion of the total variance of the time series due 
to the corresponding principal components for the optovestibular nystagmus. Panel 
e: First principal component for the vestibular nystagmus. Panel f. second prin- 
cipal component for the vestibular nystagmus time series. Panel g: third principal 
component for the vestibular nystagmus time series. Panel h: fourth principal com- 
ponent for the vestibular nystagmus time series. Panel i: fifth principal component 
for the vestibular nystagmus time series. 
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